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I.  INTRODUCTION 


For  a  mm  bar  of  alactro-optical  system  applications  (e.g.,  laser  communi¬ 
cations)  it  is  dasirabla  to  have  available  quantitative  estimates  of  various 
statistical  quantities  that  are  associated  with  atmospheric  turbulence-induced 
irradlance  fluctuations  (  i  •  e. ,  scintillation).  To  this  end,  we  have  calcu¬ 
lated  the  magnitude,  temporal  power  spectrum,  and  fade/aurge  statistics  for 
the  optical  scintillation  of  a  monochromatic  ground-based  laser  source  operat¬ 
ing  in  the  1  to  10  isi  wavelength  region,  as  observed  by  a  small  receiving 
aperture  at  an  exo-atmospherlc  location. 

As  indicated  below,  scintillation  effects  are  primarily  the  result  of 

high-altitude  dear  air  turbulence.  For  ground- to-space  applications  one  must 

employ  an  appropriate  index  structure  constant  profile  model  in  order  to 

obtain  numerical  results.  Although  a  number  of  C2  models  exist  in  the  liter*- 

Q 

ture,  we  will  apply  one  developed  by  Bufnagel.^^  In  addition  to  being  mathe¬ 
matically  convenient,  this  model  provides  fairly  good  agreement  between  the 
observed  and  predicted  values  of  the  log-amplitude  variance  of  irradlance  from 
stellar  sources.  Furthermore,  as  discussed  below,  Bufnagel's  model  agrees 

fairly  well  with  measured  values  of  C2  profiles  (e.g.,  see  Table  VII  -  IX  of 

n 

Ref.  2).  In  any  event,  although  the  numerical  results  presented  here  are 
baaed  on  the  Bufnagel  model,  the  analytic  expressions  presented  below  are 
valid  for  an  arbitrary  profile.  This  turbulence  model  contains  one  in¬ 
dependent  parsmeter,  the  root-mean-square  wind  speed  V  between  altitudes  of  5 
and  20  km  above  mean  sea  level,  l.e. , 

▼  -  [jy  /20  v2(h)dh)]1/2. 
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The  quantity  V  can  ba  obtaiaad,  for  example,  from  rawlndaonda  data  and  la 


typically  in  tha  20  -  35  a/s  range.  In  Ref.  1.  Hufnagal  coanants  that,  for 
the  State  of  Maryland,  V  la  noraally  distributed  about  27  n/a  with  a  standard 
deviation  of  9  a/s.  Furthermore,  for  gTound-t ©-space  propagation  conditions 
the  characteristic  irradlance  correlation  length  at  the  location  of  the 
receiver  is  typically  auch  larger  than  tha  dlaaater  of  the  receiving  aper¬ 
ture.  As  such,  we  can  aasuae  that  tha  receiver  is  a  "point-receiver"  and  thus 
no  aperture  averaging  effects  are  expected. 


II.  SPATIAL  STATISTICS 

For  all  cases  of  interest  here,  the  variance  of  Irradlance  when  normal¬ 
ised  by  the  aean  irradlance  level  Is  auch  lass  than  unity.  Therefore,  the  log 
irradlance  1  is  noraally  distributed.  It  is  operationally  useful  to  define 
the  log-lrradlance  1  by 

i(r,t)  -  — j - i  e\  (1) 


where  I<£,t)  la  the  irradlance  at  position  vector  _£_  and  time  t,  and  Ig  Is  Its 
aean  value.  If  the  variance  of  1  is  snail  compared  to  unity,  the  probability 
distribution  for  i  Is  normal 


Pt(i) 


,  u  -  i  r 

«*p[ - 5 — ], 

2 «  * 


where  i  la  the  aean  value  of  1,  and 
■ 


(2) 


o  2  -  «l  -  1  )2>  -  <[ln[l/Ia)]2>. 


(3) 


Also,  it  can  be  shown  that  A  ■  <A>  »  -( l/2)o .  .  Then,  the  probability 

a  * 

2 

distribution  for  log  irradiance  is  characterised  by  one  quantity,  er£  ,  the 
variance  of  log  irradiance: 


PA(i>  - 


<2™A  > 


f  [A  ♦  d/2)o/] 

exp  i - * - ] . 


It  follows  that  the  distribution  function  for  PA(i)  is  log  normal  and  is  given 


P.U)  - \  77 7  t-  e*p{ - ^  [In  i  +  i  o  2]  }. 

1  (2»o/)1/2  1  2c/  2  * 


While  o/,  the  variance  of  log-lrradiance ,  appears  in  the  literature  as  the 
most  common  measure  of  the  magnitude  of  scintillation,  it  can  be  simply 
related  to  the  normalized  irradiance  variance  as 


2  <(I  ’  Im)>  r  2i  2  2 

c/  s - J-S - -  exp[o/]  -  1  -  c/,  for  o/«l. 


2 

To  proceed,  estimates  of  o^  are  needed .  From  Ref.  3  it  can  be  shown  that 


.  2.2*  k7'6(..c9)u'6  r  c  1  <.>  h5/6dh, 

x>  n 

o 


where  0  is  the  zenith  angle  (for  the  case  at  hand  9  <  60°),  k  is  the  optical 

2 

wave  number  (-  2*/A),  H  is  the  altitude  of  the  receiver,  and  C  (h)  Is  the 

n 

index  structure  constant  profile,  i.e.,  the  strength  of  turbulence  as  a  func¬ 
tion  of  height  above  ground.  Because  of  the  presence  of  the  weighting 

2 

function  in  the  Integrand  of  Eq.  (7),  o^  is  very  insensitive  to  values 
2 

of  C  close  to  the  ground.  Rather,  it  is  the  high-altitude  turbulence  in  the 


tropopauae  at  an  altitude  of  about  9  to  12  ka  that  Is  aost  effective  in  pt  >- 
duclng  optical  scintillation  from  a  source  on  the  ground. Furthermore,  for 
receiver  altitudes  H  well  above  the  turbulent  ataosphere  (i.e.,  H  »  20  ka), 
the  value  of  the  integral  appearing  in  Eq.  (7)  is  independent  of  H,  from  which 
it  follows  that  the  resulting  log-lrradiance  variance  does  not  depend  on  the 
altitude  of  the  receiver. 

On  the  basis  of  various  eaplrlcal  scintillation  data,  Hufnagel  postulated 
that  for  optical  and  infrared  wavelengths,  the  index  of  refraction  structure 
constant  for  altitudes  greater  than  3  ka  above  aean  sea  level  is(l) 

C2(h)  -  2.72  x  10“16[3  V2  Ojjr10  ®“h  ♦  exp  (-h/1.5)]  (m“2/3),  (8) 

where  V2  is  the  Man  value  of  the  wind  speed  squared  in  meters  squared  per 

seconds  squared  in  the  5  to  20  ka  range,  and  h  is  the  altitude  in  kilo** 

asters.  Recent  measurements  have  confined  the  approximate  validity  of  this 

aodel.  Barlettl  et  al.*  measured  C2  at  three  different  sites  in  Europe  and 

found  that  for  h  >  4  ka  their  results  were  nearly  independent  of  site  and 

agreed  fairly  well  with  the  Hufnagel  aodel.  Vernln  et  al.5  measured  C2  at 

n 

numerous  sites  in  Italy,  and  again  their  results  did  not  deviate  greatly  from 
the  corresponding  predicted  values  of  the  aodel.  These  measurements  also 
indicated  that  there  was  very  little  seasonal  variation  in  C‘,  the  dominant 
quantity  in  all  cases  being  the  value  of  the  root -assn-square  wind  speed  V. 

Substitution  of  Eq.  (8)  into  Eq.  (7)  and  performance  of  the  integration 
yields 

«  k/.41  x  10“2  C^)2  +  4.45  *  10”3]  X”7/6(sec9)11/6,  (9) 


where  X  Is  in  units  of  l*.  Figure  1  is  a  plot,  of  Eq.  (9)  at  X  *  2.91  )ta  as  a 

function  of  V  for  various  values  of  6.  For  example,  at  0  -  60°  and  V  »  30 

m/s,  o^  *  9.8  *  10  .  In  fact,  for  all  cases  of  practical  interest,  (V  <  50 

o 

m/s)  o^  ,  as  measured  above  the  earth's  atmosphere  and  looking  within  zenith 

2 

angles  of  less  than  60 is  always  much  less  than  unity. 

In  Table  1  the  log-lrradiance  variance  is  given  for  various  Infrared 
wavelengths  of  Interest  for  two  representative  root-mean-square  wind  speed 
values.  The  two  values  of  V  used  represent,  for  the  State  of  Maryland,  the 
swan  and  one  standard  deviation  above  the  mean  value.  As  such,  these  cases 
should  represent  nominal  and  fairly  strong  turbulence  conditions.  Although 
larger  values  of  V  may  be  appropriate,  they  are  relatively  unlikely.  In  any 

event,  if  larger  values  of  V  are  deemed  operationally  applicable,  the  corres¬ 

ponding  numerical  results  can  be  obtained  in  a  straightforward  manner  from  the 
analytic  results  presented  below. 

Another  quantity  of  Interest  is  the  covariance  of  normalized  irradiance 
Cj,  which  is  defined  as 


Ci(-r-2)  "  <[1^Ei»t J  “  l-ia ] [ 1  [ E2 * c ^  "  1ml>»  (W) 

Here  r^and  ^  ere  two  points  in  a  plane  perpendicular  to  the  line  of  sight  at 
the  receiver  altitude.  Assuming  the  covariance  of  irradiance  Is  stationary  in 
the  wide  sense,  it  follows  that  C  (jj.r^)  *  Ci^  where  r  -  r^  r It  can 
be  shown  that  ig  related  to  the  corresponding  covariance  of  log-lrradiance 

C*a* 

C  (r)  -  exp[C  (r)]  -1  -  C.(r),  for  a  2  «  1, 


(11) 


Table  1.  Log-irradlance  variance  at  several  wavelengths  for 
various  values  of  the  root-aaan-square  wind  speed  and  zenith  angles . 


B 


where 


Ct(r)  -  <[i(r1)  -  laIti(r2)  -  ij>.  (12) 

Hoc*  that  0^(0)  ■  Hence,  w«  writ*  C^Cr)  -  o^^b^Cr),  where  b£(r)  Is  tha 

notullnd  irradianca  correlation  coefficient. 

It  can  b«  shown  that^*»^ 

/  C2  (h)  dh  /  J0[Kr(h/H)]  sin^lAi/lk]*  "8/3dK 

bA(r)  -  -2 - g - 2 - - -  ,  (13) 

/  C2  (h)  dh  /  sin2  [k2!^]  K~8/3dK 
o  o 

vhsrs  Jq  la  Cha  zaro-order-based  function  of  tha  first  kind,  tfa  do  not  dwell 
upon  this  quantity  at  any  langth  azcapt  to  nota  that  tha  charactaristlc  scala 
length  r^,  daflnad  as  that  value  of  r  such  that  b^r^)  ■  0J  is,  for  tha  case 
at  hand,  of  tha  ordar  of 


r,  -  (£<W« 

where  hQ  ■  10  -  12  ka,  and  dh  -  3  -  5  ka.  Tha  precise  values  of  1^  and  Ah  de¬ 
pend  weakly  on  V  but  are  always  in  tha  range  given  above  for  all  values  of  V 
of  practical  concern*  Hence*  in  the  Infrared  and  for  H  »  200 

ka,  r£  »  1  a,  which  is  auch  larger  than  the  receiver  aperture  assuaed  here. 
Thus,  we  can  assuaa  that  the  receiver  is  a  "point"  receiver,  and  no  aperture 
averaging  effects  are  expected* 

III.  TEMPORAL  STATISTICS 

The  teaporal  changes  in  intensity  are  obtained  froa  the  assumption  that 
tha  wind  sweeps  different  turbulent  eddies  past  the  line  of  sight  as  a  func- 


tion  of  time.  This  is  th«  "frozen  in'*  nodal  of  turbulence  and  la  justified 
for  the  present  calculation.  Of  primary  interest  is  the  temporal  power  spec¬ 
trum  of  irradiance  or.  equivalently,  log-lrradiance .  As  long  as  the  log- 
lrradience  variance  is  much  less  than  unity,  as  it  is  in  the  present  case,  the 
irradiance  power  spectrum  is  equal  to  the  log-irradiance  power  spectrum.  For 
convenience,  we  present  our  results  in  terms  of  the  log-irradiance  power 
spectrum.  The  temporal  power  spectrum  P(f)  is  defined  as  the  Fourier  trans¬ 
form  of  the  covariance  function,  and  for  the  present  case,  it  can  be  shown 
that*3* 


7/,  H  C2(h)h4/3 

P(f)  «  8.27  (sec9)7/V/3  /  -H__- dh 

o  n 

2  .  f2  i-H/6  ._2r_2  .  .2, .2 


(15) 


*  /  [*'  + 
o 


*?h) 


sin  [x*  +  fVf“(h)]dx 


where  vQ(h)  is  the  wind  speed  normal  to  the  direction  of  propagation  at  alti¬ 
tude  h,  and 

*.<w  ’  <16> 

Equation  (15)  has  been  integrated  numerically  for  various  cases  of  vQ(h) , 
obtained  from  simulated  wind  speed  data,  and  the  results  are  plotted  in  Figure 
2  for  various  values  of  V,  the  root-mean-square  wind  speed,  between  5  and 
20  km.  Examination  of  Figure  2  reveals  that  the  bandwidth  of  the  fluctuations 
are  of  the  order  of  several  hundred  hertz  for  all  cases  of  Interest. 


IV.  FADE  ADD  SUSGS  STATISTICS 


S«xt,  the  results  for  various  probabilities  of  fad* a  and  surges  are  pre- 
s*nc*d .  In  vfaac  follows,  all  lrradlancs  values  are  assumed  to  be  normalized 
to  the  mean  value  IB.  It  Is  customary  to  refer  to  fades  or  surges  In 
decibels.  That  Is,  If  the  value  of  lrradlance  is  I,  where  I  is  greater  (less) 
than  1^,  then  we  are  dealing  with  a  surge  (fade)  of  magnitude  S  dB,  (F  dB) , 
where 

S  -  10  log(I/I  )  (17) 

11 

and 

F  -  -  10|log(I/IB)|  (18) 

The  relationship  between  a  given  fade  or  surge  level  and  the  corresponding 
level  of  log-lrradlance  Is  obtained  from  Eqs.  (1),  (17),  and  (18)  as 

-  [Jtn(10)/10]S  -  0.23  S  (Surges)  (19a) 

and 

Ip  -  -[In(10)/10JF  -  -0.23  |?|  (Fades)  (19b) 

We  also  present  results  for  the  following:  (1)  the  fraction  of  tine  that 
the  lrradlance  is  above  or  below  a  given  value,  (2)  the  aean  number  of  tines 
per  second  that  the  lrradlance  crosses  a  given  value,  and  (3)  the  aean  dura¬ 
tion  of  fade  or  a  surge  at  a  given  level  I.  Again,  we  let  i(t)  denote  the 
instantaneous  value  of  lrradlance,  normalized  to  its  aean. 
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A.  Fraction  of  Tine  of  Fades  and  Surges 


We  obtain  the  fraction  of  tine  that  i(t)  >  1  and  i(t)  <  i  ,  respective - 

o  o 

ly.  Since  we  are  dealing  with  an  ergodic  process,  where  time  averages  are 

equal  to  enseable  averages,  the  fraction  of  time  of  a  surge  where  1  >  1  [de- 

o 

noted  by  Frac(i  >  i  )]  is 

o 

Frac(i  >  i  )  -  P,(i  >  i  )  -  P,(i  >  *  ),  (20) 

ox  o  *  o 

where  P^  *nd  are  the  cunulatlve  probabilities  for  irradlance  and  log-lrrad- 
lance,  respectively,  and  1q  ■  exp(ifl).  Similarly,  the  corresponding  fraction 
of  tine  of  a  fade  Frac(i  <  iQ)  is 


Frac(i  <  i  )  -  P4(i  <  i  )  -  Pt(i  <  l  ). 

OX  o  *  o 


(21) 


Since  i  is  normally  distributed,  we  obtain 


F~c(i  <  V  -  73TT7I  J. 


(i  +  1/2  oh2 
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.  wa/ 

/  exp[ - 5 - ]  di 
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where  erf  denotes  the  error  function. 


-  0.23  F 
Po  ° 

and  F0  1*  the  fade  level  in  decibels.  Sinllarly,  we  obtain 


(23) 
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and  SQ  la  the  surge  level  in  decibels •  Figures  3  and  4  are  plots  of  Eqs .  (22) 

2 

and  (24)  as  a  function  of  o^  for  various  fade  and  surge  levels. 


B.  Mean  Frequency  of  Level  Crossings 

Another  quantity  of  Interest  Is  the  frequency  of  fades  and  surges.  This 
Is  obviously  Identical  to  the  frequency  of  positive  or  negative  crossings  of 
the  level  JtQ  (l.e.,  1Q  -  Iq/I*),  positive  and  negative  referring  to  the  slope 
of  Kt)  at  the  tins  of  crossing.  We  denote  by  v"(to)  the  nean  frequency  of 
negative  crossings.  As  a  result  of  synnatry,  the  nean  frequency  of  positive 
crossings  at  level  iQ  Is  also  equal  to  v  dQ).  Since  1  is  normally  distri- 
bated  with  a  nean  equal  to  -  (1/2)0^  and  variance  o^  ,  It  can  be  shown''0' 


[*„  +  U/»)«.2]2, 

v  <V  ■  - ^2 - 1 


where 


v  -  [/  f2P(f)df//  P(f)df]1/2 


and  P(f)  Is  the  log-irradiance  power  spectrum.  The  quantity  vq  Is  referred  to 
In  the  literature  as  the  quasi-frequency  We  numerically  Integrated  Eq. 


V.v.  -V  \  '  -V-'  v  ov  . 


2 

(27)  and  obtain  that  v  la  a  relatively  Insensitive  function  of  0.  for  values 

o  * 

2 

of  0g  of  practical  concern.  We  obtain  that  at  X  »  2.91  um ,  is  of  the 

order  530  to  550  Ha  over  the  range  of  Interest.  Hence,  for  a  fade  of  FQ  dB, 
we  obtain 

2  2 

,  1-0.23  F  ♦  (1/2)0.*] 

v'(F  )  •  Vo  [ - - 5 -  ]  (28) 
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Similarly,  for  *  surge  of  SQ  dB,  we  obtain  that  the  mean  frequency  of  positive 
crossings  is 


v*(So) 


vq  exp[- 


10.23  SQ  ♦  0/2)0^ 

20.2 


(29) 


2 

For  eugtle,  Eq.  (28)  is  plotted  in  Fig.  5  as  a  function  of  o^  for  various 

2 

fade  values.  Note  the  is,  strictly  speaking,  a  function  of  o^  .  However, 

2 

as  previously  wentloned,  it  is  a  relatively  insensitive  function  of  o^  ,  and 

to  obtain  a  worst-case  estimate,  we  use  the  value  of  v  ■  550  Hz 

o 

at  X  ■  2.91  ua.  Values  of  v  at  othsr  wavelengths  of  interest  can  be  obtained 

o 

-1/2 

directly  by  noting  that  is  proportional  to  X 

In  Figs.  6  and  7,  we  plot  the  mean  level  crossing  frequency  as  a  function 
of  fade  level  for  various  wavelengths  of  interest  in  the  infrared  for  zenith 
angles  of  0  and  60*,  respectively.  In  Fig.  8,  we  plot  the  level  crossing 
frequency  as  a  function  of  fade  level  for  X  -  1.32  ua  and  various  values  of 
the  root -me  an -square  wind  speed  and  zenith  angles.  An  examination  of  Fig.  7, 
for  example,  reveals  that  the  aean  level  crossing  frequency  associated  with  a 
3  dB  fade  is  about  70-75  Hz  at  3.8  ua  and  a  60°  zenith  angle. 
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C.  Hun  Duration  of  Fades  and  Surges 


The  aean  duration  of  a  fade  at  level  l  is 

o 


T'U#> 
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v-  U0) 


to  J3  F„  *  (1/2)0,  ] 

«p - 2-j - — 
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}]  (30) 
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and  is  plotted  in  Fig.  9  as  a  function  of  o^  •  However,  as  previously  men¬ 
tioned,  it  is  relatively  insensitive  function,  and  to  obtain  a  worst-case 
estimate,  we  use  the  value  vq  -  330  Hz.  Siailarly,  the  aean  duration  of  a 
surge  at  level  iQ  is 


T+U0) 


>  V 


(31) 
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V.  DISCUSSION 

Statistical  estiaates  of  the  coastuni cation  paraaeters  fade,  surge,  fade 
duration,  surge  duration,  frequency  of  surge,  and  frequency  of  fade  were 
derived  for  Infrared  laser  ground-to-space  communication  systeas  operating  in 
the  1  to  10  ia  wavelength  region.  The  equations  are  given  in  a  fora  applic¬ 
able  for  aost  ataospherlc  conditions  of  practical  Interest,  e.g.,  a  root-aaan- 
squara  wind  speed  of  36  a/s,  with  X  «  2.91  ua  and  a  llne-of-sight  zenith  angle 
of  60* •  For  these  conditions,  we  obtain  froa  Fig.  1  (or  Table  1) 
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Fig.  9.  Mean  duration  of  3  and  6  dB  fades  as  a  function  of 


noraallzed  log  irradiance  variance. 


Chat  *  0.14.  Once  a £  is  determined  for  the  appropriate  root -mean-square 
wind  speed,  we  are  then  In  a  position  to  obtain  the  fade  and  surge  statistics 
directly.  We  find  from  Fig.  3  that  the  fraction  of  time  the  irradlance  is 
3  dB  below  the  mean  Intensity  is  *  0.05.  Correspondingly,  the  mean  frequency 
and  duration  of  a  3  dB  fade  is,  ,  from  Figs.  5  and  9,  150  Hz  and  7  x  io~4  s, 
respectively.  The  statistics  of  other  cases  can  be  calculated  with  equal 
rapidity  and,  when  used  in  conjunction  with  other  Information,  enable  one  to 
estimate  the  availability  of  a  given  communication  channel. 


LABORATORY  OPERATIONS 


The  Laboratory  Operations  of  The  Aerospace  Corporation  la  conducting  exper¬ 
imental  and  theoretical  Investigations  necessary  for  the  evaluation  and  applica¬ 
tion  of  scientific  advances  to  naw  military  apace  systems.  Versatility  and 
flexibility  have  been  developed  to  a  high  degrae  by  the  laboratory  personnel  in 
dealing  with  the  many  problams  ancountered  In  the  nation's  rapidly  developing 
space  systems.  Expertise  In  the  latest  scientific  developments  is  vital  to  the 
accomplishment  of  tasks  related  to  these  problems.  The  laboratories  that  con¬ 
tribute  to  this  research  are: 

Aerophyslcs  Laboratory:  Launch  vehicle  and  reentry  aerodynamics  and  heat 
transfer,  propulsion  chemistry  and  fluid  mechanics,  structural  mechanics,  flight 
dynamics;  hlgh-tamperature  thermomechanica,  gas  kinetics  and  radiation;  research 
In  environmental  chemistry  end  contamination;  ew  and  pulsed  chemical  lamer 
development  Including  chemical  kinetics,  spectroscopy,  optical  resonators  and 
beam  pointing,  atmospheric  propagation,  laser  effects  and  countermeasures. 

Chemistry  and  Physics  Laboratory:  Atmospheric  chemical  reactions,  atmo- 
spherlc  optica,  light  scattering,  state-specific  chemical  reactions  and  radia¬ 
tion  transport  In  rocket  plumes,  applied  laser  spectroscopy,  laser  chemistry, 
battery  electrochemistry,  space  vacuum  and  radiation  effects  on  materials,  lu¬ 
brication  and  surface  phenomena,  thermionic  amission,  photosensitive  materials 
and  detectors,  atomic  frequency  standards,  and  bloenvlronmental  research  and 
monitoring. 

Electronics  Research  Laboratory:  Microelectronics,  GaAe  low-noise  and 
power  devices,  semiconductor  lasers,  electromagnetic  and  optical  propagation 
phenomena,  quantum  electronlca,  laser  communications,  Ildar,  and  electro-optics; 
communication  sciences,  applied  electronics,  semiconductor  crystal  and  device 
physics,  radiometric  Imaging;  millimeter-wave  and  microwave  technology. 

Information  Sciences  Research  Office:  Program  verification,  program  trans¬ 
lating  performance- sensl t lve  system  Teslgn,  distributed  architectures  for 
apaceborne  computers,  fault-tolerant  computer  systems,  artificial  Intelligence, 
and  microelectronics  applications. 

Materials  Sciences  Laboratory;  Development  of  new  materials:  metal  matrix 
composites,  polymers,  snd  new  lores  of  carbon;  component  failure  analysis  and 
reliability;  fracture  mechanics  and  stress  corrosion;  evaluation  of  materials  in 
space  environment;  materials  performance  In  space  transportation  systems;  anal¬ 
ysis  of  systems  vulnerability  and  survivability  In  enemy-induced  environments. 

Space  Sciences  Laboratory:  Atmospheric  and  Ionospheric  physics,  radiation 
from  the  atmosphere,  density  and  composition  of  the  upper  atmosphere,  aurorae 
and  alrglow;  magnetospherlc  physics,  cosmic  rays,  generation  and  propagation  of 
plasma  waves  la  the  magnetosphere;  solar  phyeles,  Infrared  astronomy;  the 
effects  of  nuclear  explosions,  aaynetlc  storms,  and  solar  activity  on  the 
earth's  atmosphere,  Ionosphere,  end  magnetosphere;  the  effects  of  optical, 
electromagnetic,  and  particulate  radiations  In  space  on  space  systems. 


